Stepwise introduction of stabilizing mutations reveals nonlinear additive effects in de novo TIM barrels.

Over the past decades, the TIM-barrel fold has served as a model system for the exploration of how changes in protein sequences affect their structural, stability, and functional characteristics, and moreover, how this information can be leveraged to design proteins from the ground up. After numerous attempts to design de novo proteins with this specific fold, sTIM11 was the first validated de novo design of an idealized four-fold symmetric TIM barrel. Subsequent efforts to enhance the stability of this initial design resulted in the development of DeNovoTIMs, a family of de novo TIM barrels with various stabilizing mutations. In this study, we present an investigation into the biophysical and thermodynamic effects upon introducing a varying number of stabilizing mutations per quarter along the sequence of a four-fold symmetric TIM barrel. We compared the base design DeNovoTIM0 without any stabilizing mutations with variants containing mutations in one, two, three, and all four quarters-designated TIM1q, TIM2q, TIM3q, and DeNovoTIM6, respectively. This analysis revealed a stepwise and nonlinear change in the thermodynamic properties that correlated with the number of mutated quarters, suggesting positive nonadditive effects. To shed light on the significance of the location of stabilized quarters, we engineered two variants of TIM2q which contain the same number of mutations but positioned in different quarter locations. Characterization of these TIM2q variants revealed that the mutations exhibit varying effects on the overall protein stability, contingent upon the specific region in which they are introduced. These findings emphasize that the amount and location of stabilized interfaces among the four quarters play a crucial role in shaping the conformational stability of these four-fold symmetric TIM barrels. Analysis of de novo proteins, as described in this study, enhances our understanding of how sequence variations can finely modulate stability in both naturally occurring and computationally designed proteins.

Exploring the Protein Sequence Space with Global Generative Models.

Recent advancements in specialized large-scale architectures for training images and language have profoundly impacted the field of computer vision and natural language processing (NLP). Language models, such as the recent ChatGPT and GPT-4, have demonstrated exceptional capabilities in processing, translating, and generating human language. These breakthroughs have also been reflected in protein research, leading to the rapid development of numerous new methods in a short time, with unprecedented performance. Several of these models have been developed with the goal of generating sequences in novel regions of the protein space. In this work, we provide an overview of the use of protein generative models, reviewing (1) language models for the design of novel artificial proteins, (2) works that use non-transformer architectures, and (3) applications in directed evolution approaches.

Retracing the evolution of a modern periplasmic binding protein.

Investigating the evolution of structural features in modern multidomain proteins helps to understand their immense diversity and functional versatility. The class of periplasmic binding proteins (PBPs) offers an opportunity to interrogate one of the main processes driving diversification: the duplication and fusion of protein sequences to generate new architectures. The symmetry of their two-lobed topology, their mechanism of binding, and the organization of their operon structure led to the hypothesis that PBPs arose through a duplication and fusion event of a single common ancestor. To investigate this claim, we set out to reverse the evolutionary process and recreate the structural equivalent of a single-lobed progenitor using ribose-binding protein (RBP) as our model. We found that this modern PBP can be deconstructed into its lobes, producing two proteins that represent possible progenitor halves. The isolated halves of RBP are well folded and monomeric proteins, albeit with a lower thermostability, and do not retain the original binding function. However, the two entities readily form a heterodimer in vitro and in-cell. The x-ray structure of the heterodimer closely resembles the parental protein. Moreover, the binding function is fully regained upon formation of the heterodimer with a ligand affinity similar to that observed in the modern RBP. This highlights how a duplication event could have given rise to a stable and functional PBP-like fold and provides insights into how more complex functional structures can evolve from simpler molecular components.

Thermodynamic and kinetic analysis of the LAO binding protein and its isolated domains reveal non-additivity in stability, folding and function.

Substrate-binding proteins (SBPs) are used by organisms from the three domains of life for transport and signalling. SBPs are composed of two domains that collectively trap ligands with high affinity and selectivity. To explore the role of the domains and the integrity of the hinge region between them in the function and conformation of SBPs, here, we describe the ligand binding, conformational stability and folding kinetics of the Lysine Arginine Ornithine (LAO) binding protein from Salmonella thiphimurium and constructs corresponding to its two independent domains. LAO is a class II SBP formed by a continuous and a discontinuous domain. Contrary to the expected behaviour based on their connectivity, the discontinuous domain shows a stable native-like structure that binds l-arginine with moderate affinity, whereas the continuous domain is barely stable and shows no detectable ligand binding. Regarding folding kinetics, studies of the entire protein revealed the presence of at least two intermediates. While the unfolding and refolding of the continuous domain exhibited only a single intermediate and simpler and faster kinetics than LAO, the folding mechanism of the discontinuous domain was complex and involved multiple intermediates. These findings suggest that in the complete protein the continuous domain nucleates folding and that its presence funnels the folding of the discontinuous domain avoiding nonproductive interactions. The strong dependence of the function, stability and folding pathway of the lobes on their covalent association is most likely the result of the coevolution of both domains as a single unit.

Structures of permuted halves of a modern ribose-binding protein.

Periplasmic binding proteins (PBPs) are a class of proteins that participate in the cellular transport of various ligands. They have been used as model systems to study mechanisms in protein evolution, such as duplication, recombination and domain swapping. It has been suggested that PBPs evolved from precursors half their size. Here, the crystal structures of two permuted halves of a modern ribose-binding protein (RBP) from Thermotoga maritima are reported. The overexpressed proteins are well folded and show a monomer-dimer equilibrium in solution. Their crystal structures show partially noncanonical PBP-like fold type I conformations with structural deviations from modern RBPs. One of the half variants forms a dimer via segment swapping, suggesting a high degree of malleability. The structural findings on these permuted halves support the evolutionary hypothesis that PBPs arose via a duplication event of a flavodoxin-like protein and further support a domain-swapping step that might have occurred during the evolution of the PBP-like fold, a process that is necessary to generate the characteristic motion of PBPs essential to perform their functions.

A newly introduced salt bridge cluster improves structural and biophysical properties of de novo TIM barrels.

Protein stability can be fine-tuned by modifying different structural features such as hydrogen-bond networks, salt bridges, hydrophobic cores, or disulfide bridges. Among these, stabilization by salt bridges is a major challenge in protein design and engineering since their stabilizing effects show a high dependence on the structural environment in the protein, and therefore are difficult to predict and model. In this work, we explore the effects on structure and stability of an introduced salt bridge cluster in the context of three different de novo TIM barrels. The salt bridge variants exhibit similar thermostability in comparison with their parental designs but important differences in the conformational stability at 25°C can be observed such as a highly stabilizing effect for two of the proteins but a destabilizing effect to the third. Analysis of the formed geometries of the salt bridge cluster in the crystal structures show either highly ordered salt bridge clusters or only single salt bridges. Rosetta modeling of the salt bridge clusters results in a good prediction of the tendency on stability changes but not the geometries observed in the three-dimensional structures. The results show that despite the similarities in protein fold, the salt bridge clusters differently influence the structural and stability properties of the de novo TIM barrel variants depending on the structural background where they are introduced.

The Stability Landscape of de novo TIM Barrels Explored by a Modular Design Approach.

The ability to design stable proteins with custom-made functions is a major goal in biochemistry with practical relevance for our environment and society. Understanding and manipulating protein stability provide crucial information on the molecular determinants that modulate structure and stability, and expand the applications of de novo proteins. Since the (ß/⍺)8-barrel or TIM-barrel fold is one of the most common functional scaffolds, in this work we designed a collection of stable de novo TIM barrels (DeNovoTIMs), using a computational fixed-backbone and modular approach based on improved hydrophobic packing of sTIM11, the first validated de novo TIM barrel, and subjected them to a thorough folding analysis. DeNovoTIMs navigate a region of the stability landscape previously uncharted by natural TIM barrels, with variations spanning 60 degrees in melting temperature and 22 kcal per mol in conformational stability throughout the designs. Significant non-additive or epistatic effects were observed when stabilizing mutations from different regions of the barrel were combined. The molecular basis of epistasis in DeNovoTIMs appears to be related to the extension of the hydrophobic cores. This study is an important step towards the fine-tuned modulation of protein stability by design.

Evolution, folding, and design of TIM barrels and related proteins.

Proteins are chief actors in life that perform a myriad of exquisite functions. This diversity has been enabled through the evolution and diversification of protein folds. Analysis of sequences and structures strongly suggest that numerous protein pieces have been reused as building blocks and propagated to many modern folds. This information can be traced to understand how the protein world has diversified. In this review, we discuss the latest advances in the analysis of protein evolutionary units, and we use as a model system one of the most abundant and versatile topologies, the TIM-barrel fold, to highlight the existing common principles that interconnect protein evolution, structure, folding, function, and design.

The interplay of protein-ligand and water-mediated interactions shape affinity and selectivity in the LAO binding protein.

The study of binding thermodynamics is essential to understand how affinity and selectivity are acquired in molecular complexes. Periplasmic binding proteins (PBPs) are macromolecules of biotechnological interest that bind a broad number of ligands and have been used to design biosensors. The lysine-arginine-ornithine binding protein (LAO) is a PBP of 238 residues that binds the basic amino acids l-arginine and l-histidine with nm and µm affinity, respectively. It has been shown that the affinity difference for arginine and histidine binding is caused by enthalpy, this correlates with the higher number of protein-ligand contacts formed with arginine. In order to elucidate the structural bases that determine binding affinity and selectivity in LAO, the contribution of protein-ligand contacts to binding energetics was assessed. To this end, an alanine scanning of the LAO-binding site residues was performed and arginine and histidine binding were characterized by isothermal titration calorimetry and X-ray crystallography. Although unexpected enthalpy and entropy changes were observed in some mutants, thermodynamic data correlated with structural information, especially, the binding heat capacity change. We found that selectivity is conferred by several residues rather than exclusive arginine-protein interactions. Furthermore, crystallographic structures revealed that protein-ligand contributions to binding thermodynamics are highly influenced by the solvent. Finally, we found a similar backbone conformation in all the closed structures obtained, but different structures in the open state, suggesting that the binding site residues of LAO play an important role in stabilizing not only the holo conformation, but also the apo state. DATABASE: Structural data are available in the Protein Data Bank database under the accession numbers 6MLE, 6MLN, 6MLG, 6MKX, 6MLI, 6MLA, 6MKU, 6MKW, 6ML0, 6MLD, 6MLV, 6MLO, 6MLP, 6ML9, 6MLJ.

Voltage vs. Ligand I: Structural basis of the intrinsic flexibility of S3 segment and its significance in ion channel activation.

We systematically predict the internal flexibility of the S3 segment, one of the most mobile elements in the voltage-sensor domain. By analyzing the primary amino acid sequences of V-sensor containing proteins, including Hv1, TPC channels and the voltage-sensing phosphatases, we established correlations between the local flexibility and modes of activation for different members of the VGIC superfamily. Taking advantage of the structural information available, we also assessed structural aspects to understand the role played by the flexibility of S3 during the gating of the pore. We found that S3 flexibility is mainly determined by two specific regions: (1) a short NxxD motif in the N-half portion of the helix (S3a), and (2) a short sequence at the beginning of the so-called paddle motif where the segment has a kink that, in some cases, divide S3 into two distinct helices (S3a and S3b). A good correlation between the flexibility of S3 and the reported sensitivity to temperature and mechanical stretch was found. Thus, if the channel exhibits high sensitivity to heat or membrane stretch, local S3 flexibility is low. On the other hand, high flexibility of S3 is preferentially associated to channels showing poor heat and mechanical sensitivities. In contrast, we did not find any apparent correlation between S3 flexibility and voltage or ligand dependence. Overall, our results provide valuable insights into the dynamics of channel-gating and its modulation.

Voltage vs. Ligand II: Structural insights of the intrinsic flexibility in cyclic nucleotide-gated channels.

In the preceding article, we present a flexibility analysis of the voltage-gated ion channel (VGIC) superfamily. In this study, we describe in detail the flexibility profile of the voltage-sensor domain (VSD) and the pore domain (PD) concerning the evolution of 6TM ion channels. In particular, we highlight the role of flexibility in the emergence of CNG channels and describe a significant level of sequence similarity between the archetypical VSD and the TolQ proteins. A highly flexible S4-like segment exhibiting Lys instead Arg for these membrane proteins is reported. Sequence analysis indicates that, in addition to this S4-like segment, TolQ proteins also show similarity with specific motifs in S2 and S3 from typical V-sensors. Notably, S3 flexibility profiles from typical VSDs and S3-like in TolQ proteins are also similar. Interestingly, TolQ from early divergent prokaryotes are comparatively more flexible than those in modern counterparts or true V-sensors. Regarding the PD, we also found that 2TM K+-channels in early prokaryotes are considerably more flexible than the ones in modern microbes, and such flexibility is comparable to the one present in CNG channels. Voltage dependence is mainly exhibited in prokaryotic CNG channels whose VSD is rigid whereas the eukaryotic CNG channels are considerably more flexible and poorly V-dependent. The implication of the flexibility present in CNG channels, their sensitivity to cyclic nucleotides and the cation selectivity are discussed. Finally, we generated a structural model of the putative cyclic nucleotide-modulated ion channel, which we coined here as AqK, from the thermophilic bacteria Aquifex aeolicus, one of the earliest diverging prokaryotes known. Overall, our analysis suggests that V-sensors in CNG-like channels were essentially rigid in early prokaryotes but raises the possibility that this module was probably part of a very flexible stator protein of the bacterial flagellum motor complex.

FRET-based analysis and molecular modeling of the human GPN-loop GTPases 1 and 3 heterodimer unveils a dominant-negative protein complex.

GPN-loop GTPases 1 and 3 are required for RNA polymerase II (RNAPII) nuclear import. Gpn1 and Gpn3 display some sequence similarity, physically associate, and their protein expression levels are mutually dependent in human cells. We performed here Fluorescence Resonance Energy Transfer (FRET), molecular modeling, and cell biology experiments to understand, and eventually disrupt, human Gpn1-Gpn3 interaction in live HEK293-AD cells. Transiently expressed EYFP-Gpn1 and Gpn3-CFP generated a strong FRET signal, indicative of a very close proximity, in the cytoplasm of HEK293-AD cells. Molecular modeling of the human Gpn1-Gpn3 heterodimer based on the crystallographic structure of Npa3, the Saccharomyces cerevisiae Gpn1 ortholog, revealed that human Gpn1 and Gpn3 associate through a large interaction surface formed by internal α-helix 7, insertion 2, and the GPN-loop from each protein. In site-directed mutagenesis experiments of interface residues, we identified the W132D and M227D EYFP-Gpn1 mutants as defective to produce a FRET signal when coexpressed with Gpn3-CFP. Simultaneous but not individual expression of Gpn1 and Gpn3, with either or both proteins fused to EYFP, retained RNAPII in the cytoplasm and markedly inhibited global transcription in HEK293-AD cells. Interestingly, the W132D and M227D Gpn1 mutants that showed an impaired ability to interact with Gpn3 by FRET were also unable to delocalize RNAPII in this assay, indicating that an intact Gpn1-Gpn3 interaction is required to display the dominant-negative effect on endogenous Gpn1/Gpn3 function we described here. Altogether, our results suggest that a Gpn1-Gpn3 strong interaction is critical for their cellular function.

Structural basis for the modulation of plant cytosolic triosephosphate isomerase activity by mimicry of redox-based modifications.

Reactive oxidative species (ROS) and S-glutathionylation modulate the activity of plant cytosolic triosephosphate isomerases (cTPI). Arabidopsis thaliana cTPI (AtcTPI) is subject of redox regulation at two reactive cysteines that function as thiol switches. Here we investigate the role of these residues, AtcTPI-Cys13 and At-Cys218, by substituting them with aspartic acid that mimics the irreversible oxidation of cysteine to sulfinic acid and with amino acids that mimic thiol conjugation. Crystallographic studies show that mimicking AtcTPI-Cys13 oxidation promotes the formation of inactive monomers by reposition residue Phe75 of the neighboring subunit, into a conformation that destabilizes the dimer interface. Mutations in residue AtcTPI-Cys218 to Asp, Lys, or Tyr generate TPI variants with a decreased enzymatic activity by creating structural modifications in two loops (loop 7 and loop 6) whose integrity is necessary to assemble the active site. In contrast with mutations in residue AtcTPI-Cys13, mutations in AtcTPI-Cys218 do not alter the dimeric nature of AtcTPI. Therefore, modifications of residues AtcTPI-Cys13 and AtcTPI-Cys218 modulate AtcTPI activity by inducing the formation of inactive monomers and by altering the active site of the dimeric enzyme, respectively. The identity of residue AtcTPI-Cys218 is conserved in the majority of plant cytosolic TPIs, this conservation and its solvent-exposed localization make it the most probable target for TPI regulation upon oxidative damage by reactive oxygen species. Our data reveal the structural mechanisms by which S-glutathionylation protects AtcTPI from irreversible chemical modifications and re-routes carbon metabolism to the pentose phosphate pathway to decrease oxidative stress.

Structure and conformational stability of the triosephosphate isomerase from Zea mays. Comparison with the chemical unfolding pathways of other eukaryotic TIMs.

We studied the structure, function and thermodynamic properties for the unfolding of the Triosephosphate isomerase (TIM) from Zea mays (ZmTIM). ZmTIM shows a catalytic efficiency close to the diffusion limit. Native ZmTIM is a dimer that dissociates upon dilution into inactive and unfolded monomers. Its thermal unfolding is irreversible with a Tm of 61.6 ±â€¯1.4 °C and an activation energy of 383.4 ±â€¯11.5 kJ mol-1. The urea-induced unfolding of ZmTIM is reversible. Transitions followed by catalytic activity and spectroscopic properties are monophasic and superimposable, indicating that ZmTIM unfolds/refolds in a two-state behavior with an unfolding ΔG°(H20) = 99.8 ±â€¯5.3 kJ mol-1. This contrasts with most other studied TIMs, where folding intermediates are common. The three-dimensional structure of ZmTIM was solved at 1.8 Å. A structural comparison with other eukaryotic TIMs shows a similar number of intramolecular and intermolecular interactions. Interestingly the number of interfacial water molecules found in ZmTIM is lower than those observed in most TIMs that show folding intermediates. Although with the available data, there is no clear correlation between structural properties and the number of equilibrium intermediates in the unfolding of TIM, the identification of such structural properties should increase our understanding of folding mechanisms.

A strategy based on thermal flexibility to design triosephosphate isomerase proteins with increased or decreased kinetic stability.

Kinetic stability of proteins determines their susceptibility to irreversibly unfold in a time-dependent process, and therefore its half-life. A residue displacement analysis of temperature-induced unfolding molecular dynamics simulations was recently employed to define the thermal flexibility of proteins. This property was found to be correlated with the activation energy barrier (Eact) separating the native from the transition state in the denaturation process. The Eact was determined from the application of a two-state irreversible model to temperature unfolding experiments using differential scanning calorimetry (DSC). The contribution of each residue to the thermal flexibility of proteins is used here to propose multiple mutations in triosephosphate isomerase (TIM) from Trypanosoma brucei (TbTIM) and Trypanosoma cruzi (TcTIM), two parasites closely related by evolution. These two enzymes, taken as model systems, have practically identical structure but large differences in their kinetic stability. We constructed two functional TIM variants with more than twice and less than half the activation energy of their respective wild-type reference structures. The results show that the proposed strategy is able to identify the crucial residues for the kinetic stability in these enzymes. As it occurs with other protein properties reflecting their complex behavior, kinetic stability appears to be the consequence of an extensive network of inter-residue interactions, acting in a concerted manner. The proposed strategy to design variants can be used with other proteins, to increase or decrease their functional half-life.

Localized conformational changes trigger the pH-induced fibrillogenesis of an amyloidogenic λ light chain protein.

Solvent conditions modulate the expression of the amyloidogenic potential of proteins. In this work the effect of pH on the fibrillogenic behavior and the conformational properties of 6aJL2, a model protein of the highly amyloidogenic variable light chain λ6a gene segment, was examined. Ordered aggregates showing the ultrastructural and spectroscopic properties observed in amyloid fibrils were formed in the 2.0-8.0 pH range. At pH <3.0 a drastic decrease in lag time and an increase in fibril formation rate were found. In the 4.0-8.0 pH range there was no spectroscopic evidence for significant conformational changes in the native state. Likewise, heat capacity measurements showed no evidence for residual structure in the unfolded state. However, at pH <3.0 stability is severely decreased and the protein suffers conformational changes as detected by circular dichroism, tryptophan and ANS fluorescence, as well as by NMR spectroscopy. Molecular dynamics simulations indicate that acid-induced conformational changes involve the exposure of the loop connecting strands E and F. These results are compatible with pH-induced changes in the NMR spectra. Overall, the results indicate that the mechanism involved in the acid-induced increase in the fibrillogenic potential of 6aJL2 is profoundly different to that observed in κ light chains, and is promoted by localized conformational changes in a region of the protein that was previously not known to be involved in acid-induced light chain fibril formation. The identification of this region opens the potential for the design of specific inhibitors.

Side chain flexibility and coupling between the S4-S5 linker and the TRP domain in thermo-sensitive TRP channels: Insights from protein modeling.

The transient receptor potential (TRP) superfamily is subdivided into several subfamilies on the basis of sequence similarity, which is highly heterogeneous but shows a molecular architecture that resembles the one present in members of the Kv channel superfamily. Because of this diversity, they produce a large variety of channels with different gating and permeability properties. Elucidation of these particular features necessarily requires comparative studies based on structural and functional data. The present study aims to compilate, analyze, and determine, in a coherent way, the relationship between intrinsic side-chain flexibility and the allosteric coupling in members of the TRPV, TRPM, and TRPC families. Based on the recently determined structures of TRPV1 and TRPV2, we have generated protein models for single subunits of TRPV5, TRPM8, and TRPC5 channels. With these models, we focused our attention on the apparently crucial role of the GP dipeptide at the center of the S4-S5 linker and discussed its role in the interaction with the TRP domain, specifically with the highly-conserved Trp during this coupling. Our analysis suggests an important role of the S4-S5L flexibility in the thermosensitivity, where heat-activated channels possess rigid S4-S5 linkers, whereas cold-activated channels have flexible ones. Finally, we also present evidence of the key interaction between the conserved Trp residue of the TRP box and of several residues in the S4-S5L, importantly the central Pro. Proteins 2017; 85:630-646. © 2016 Wiley Periodicals, Inc.

Reversibility and two state behaviour in the thermal unfolding of oligomeric TIM barrel proteins.

Temperature is one of the main variables that modulate protein function and stability. Thermodynamic studies of oligomeric proteins, the dominant protein natural form, have been often hampered because irreversible aggregation and/or slow reactions are common. There are no reports on the reversible equilibrium thermal unfolding of proteins composed of (ß/α)8 barrel subunits, albeit this "TIM barrel" topology is one of the most abundant and versatile in nature. We studied the eponymous TIM barrel, triosephosphate isomerase (TIM), belonging to five species of different bacterial taxa. All of them were found to be catalytically efficient dimers. The three-dimensional structure of four enzymes was solved at high/medium resolution. Irreversibility and kinetic control were observed in the thermal unfolding of two TIMs, while for the other three the thermal unfolding was found to follow a two-state equilibrium reversible process. Shifts in the global stability curves of these three proteins are related to the organismal temperature range of optimal growth and modulated by variations in maximum stability temperature and in the enthalpy change at that temperature. Reversibility appears to correlate with the low isoelectric point, the absence of a residual structure in the unfolded state, small cavity volume in the native state, low conformational stability and a low melting temperature. Furthermore, the strong coupling between dimer dissociation and monomer unfolding may reduce aggregation and favour reversibility. It is therefore very thought-provoking to find that a common topological ensemble, such as the TIM barrel, can unfold/refold in the Anfinsen way, i.e. without the help of the cellular machinery.